Advanced video coded pictures - reduced cost computation of an intra mode decision in the frequency domain

ABSTRACT

A method of intra mode prediction uses a block of pixels and their horizontal H pos  and vertical V pos  pixel positions and adjacent horizontal and vertical pixels within an input picture frame signal as inputs for a method for selecting a lowest Sum of Absolute Transformed Differences (SATD) intra mode among intra modes among the horizontal, vertical, and steady state (DC) intra modes for use in advanced video coding algorithms, such as MPEG-4 Part 10 and H.264/AVC. Associated costs are calculated for each of the intra modes are used to select the output of the best intra mode. The method reduces the unimproved computational cost of three 2D 4×4 Hadamard transformations (which is equivalent to 24 1D 4 point transformations) to just 4 1D 4 point transformations for a significant computational improvement. As horizontal and vertical panning is frequently used in video imagery, this improvement may reduce encoder processing by 80%.

CROSS-REFERENCE TO RELATED APPLICATIONS

Not Applicable

STATEMENT REGARDING FEDERALLY SPONSORED RESEARCH OR DEVELOPMENT

Not Applicable

INCORPORATION-BY-REFERENCE OF MATERIAL SUBMITTED ON A COMPACT DISC

Not Applicable

BACKGROUND OF THE INVENTION

1. Field of the Invention

This invention pertains generally to video encoding, and more particularly to intra mode decisions within advanced video encoding (such as H.264/AVC or MPEG-4 Part 10) standards.

2. Description of Related Art

H.264/AVC, alternatively known as MPEG-4 Part 10 and by several other monikers, is representative of improved data compression algorithms. Improved data compression, however, comes at the price of greatly increased computational requirements during the encoding processing phase.

In particular, the H.264/AVC contains a very computationally expensive section where an optimal intra mode is calculated that yields the lowest sum of absolute transformed differences (or SATD) between an input 4×4 pixel block, and its compressed version. Examples of prior attempts to reduce this computational cost include A. C. Yu, G. R. Martin, and H. Park, in “A Frequency Domain Approach to Intra Mode Selection in H.264/AVC”, incorporated by reference in its entirety, U.S. patent publication US 2006/0209948 A1, incorporated herein by reference in its entirety, and U.S. patent publication US 2006/0251330 A1, incorporated herein by reference in its entirety.

BRIEF SUMMARY OF THE INVENTION

Accordingly, one aspect of the invention is a method of selecting an intra mode in a video image encoder, comprising:

(a) selecting a lowest Sum of Absolute Transformed Differences (SATD) intra mode among intra modes while using [x], {right arrow over (h)}, {right arrow over (v)}, H_(pos), V_(pos);

(b) wherein [x] is a 4×4 block of pixels within a picture and

${{\lbrack x\rbrack \equiv \begin{bmatrix} x_{0,0} & x_{0,1} & x_{0,2} & x_{0,3} \\ x_{1,0} & x_{1,1} & x_{1,2} & x_{1,3} \\ x_{2,0} & x_{2,1} & x_{2,2} & x_{2,3} \\ x_{3,0} & x_{3,1} & x_{3,2} & x_{3,3} \end{bmatrix}},{\overset{\rightarrow}{h} = \begin{bmatrix} x_{0,{- 1}} \\ x_{1,{- 1}} \\ x_{2,{- 1}} \\ x_{3,{- 1}} \end{bmatrix}},{and}}\mspace{14mu}$ ${\overset{\rightarrow}{v} = \begin{bmatrix} x_{{- 1},0} \\ x_{{- 1},1} \\ x_{{- 1},2} \\ x_{{- 1},3} \end{bmatrix}};$

(c) wherein H_(pos) is a horizontal pixel position of the 4×4 block within the image;

(d) wherein V_(pos) is a vertical pixel position of the 4×4 block within the image;

(e) wherein the lowest SATD intra mode is determined among a group comprising: (i) a horizontal intra mode; (ii) a vertical intra mode; and (iii) a steady state (DC) intra mode; and

(f) outputting the lowest SATD intra mode.

The outputting step may comprise:

(a) calculating a horizontal predictor {right arrow over (H)}≡(H₀,H₁,H₂,H₃)^(T), a vertical predictor {right arrow over (V)}≡(V₀,V₁,V₂,V₃)^(T), and a steady state (DC) predictor D;

(b) calculating a horizontal cost precursor C_(hs) and a vertical cost precursor C_(vs) using the horizontal predictor {right arrow over (H)}, the vertical predictor {right arrow over (V)}, and the steady state (DC) predictor D; and

(c) calculating a horizontal intra mode cost C_(H), a vertical intra mode cost C_(V), and a steady state (DC) intra mode cost C_(D) using the horizontal cost precursor C_(hs) and the vertical cost precursor C_(vs).

The calculating of the horizontal predictor {right arrow over (H)}, the vertical predictor {right arrow over (V)}, and the steady state (DC) predictor D may comprise:

(a) if H_(pos)≠0 and V_(pos)≠0 then:

-   -   i. setting {right arrow over (H)}≡[T₄]{right arrow over (h)}         -   wherein {right arrow over             (h)}≡(h₀,h₁,h₂,h₃)^(T)≡(x_(0,−1),x_(1,−1),x_(2,−1),x_(3,−1))^(T),     -   ii. setting {right arrow over (V)}=[T₄]{right arrow over (v)}         -   wherein {right arrow over             (V)}≡(v₀,v₁,v₂,v₃)≡(x_(−1,0),x_(−1,1),x_(−1,2),x_(−1,3))^(T),             and     -   iii. setting D=(H₀+V₀)/2;

(b) if H_(pos)=0 and V_(pos)≠0 then:

-   -   i. setting {right arrow over (H)}=(2¹⁵−1,0,0,0)^(T),     -   ii. setting {right arrow over (V)}=[T₄]{right arrow over (v)}         -   wherein {right arrow over             (v)}≡(v₀,v₁,v₂,v₃)≡(x_(−1,0),x_(−1,1),x_(−1,2),x_(−1,3))^(T),             and     -   iii. setting D=V₀;

(c) if H_(pos)≠0 and V_(pos)≠0 then:

-   -   i. setting {right arrow over (H)}=[T₄]{right arrow over (h)}         -   wherein {right arrow over             (h)}≡(h₀,h₁,h₂,h₃)^(T)≡(x_(0,−1),x_(1,−1),x_(2,−1),x_(3,−1))^(T),     -   ii. setting {right arrow over (V)}=(2¹⁵−1,0,0,0)^(T), and     -   iii. setting D=H₀;

(d) if H_(pos)=0 and V_(pos)=0 then:

-   -   i. setting {right arrow over (H)}=(2¹⁵−1,0,0,0)^(T),     -   ii. setting {right arrow over (V)}=(2¹⁵−1,0,0,0)^(T), and     -   iii. setting D=128×16.

A relatively inefficient calculation of the horizontal cost precursor C_(hs) and the vertical cost precursor C_(vs) may begin with steps comprising: (a) calculating the values X_(i,0), X_(0,i) for iε0, 1, 2, 3 using the relationship

${\lbrack X\rbrack = {{\left\lbrack T_{4} \right\rbrack \lbrack x\rbrack}\left\lbrack T_{4} \right\rbrack}^{T}},{{{where}\lbrack X\rbrack} \equiv {\begin{bmatrix} X_{0,0} & X_{0,1} & X_{0,2} & X_{0,3} \\ X_{1,0} & X_{1,1} & X_{1,2} & X_{1,3} \\ X_{2,0} & X_{2,1} & X_{2,2} & X_{2,3} \\ X_{3,0} & X_{3,1} & X_{3,2} & X_{3,3} \end{bmatrix}.}}$

This method is computationally inefficient, as only the X_(i,0), X_(0,i) for iε0, 1, 2, 3 values are required. Thus, 9 of 16 values are not required.

In another aspect of the invention, a more computationally efficient calculation of the X_(i,0), X_(0,i) for iε0, 1, 2, 3 values may comprise: providing the input [x] matrix and the [T₄] matrix; then calculating:

$\begin{bmatrix} X_{0,0} \\ X_{1,0} \\ X_{2,0} \\ X_{3,0} \end{bmatrix} = {{\left\lbrack T_{4} \right\rbrack \lbrack x\rbrack}\overset{\rightarrow}{u}}$ ${{and}\begin{bmatrix} X_{0,0} & X_{0,1} & X_{0,2} & X_{0,3} \end{bmatrix}} = {{{\overset{\rightarrow}{u}}^{T}\lbrack x\rbrack}\left\lbrack T_{4} \right\rbrack}^{T}$

where {right arrow over (u)}^(T)=[1 1 1 1].

After calculating the X_(i,0), X_(0,i) for iε0, 1, 2, 3 values, then the horizontal cost precursor

$C_{hs} = {\sum\limits_{i = 1}^{3}\; {X_{i,0}}}$

and vertical cost precursor

$C_{vs} = {\sum\limits_{j = 1}^{3}\; {X_{0,j}}}$

may be calculated.

The horizontal intra mode cost C_(H) may then be calculated by

$C_{H} = {{\sum\limits_{i = 0}^{3}\; {{H_{i} - X_{i,0}}}} + {C_{vs}.}}$

The vertical intra mode cost C_(V) may be calculated by

$C_{v} = {{\sum\limits_{j = 0}^{3}\; {{V_{j} - X_{0,j}}}} + {C_{hs}.}}$

The steady state (DC) intra mode cost C_(D) may be calculated by C_(D)=|D−X_(0,0)|+C_(hs)+C_(vs).

The calculating may further comprise: selecting the lowest SATD intra mode with a lowest associated intra mode cost among the group consisting of: the horizontal intra mode cost C_(H), the vertical intra mode cost C_(V), and the steady state (DC) intra mode cost C_(D). This calculating process is typically achieved through a series of logical operations.

Another aspect of the invention is a computer readable medium comprising a programming executable capable of performing on a computer the method described above.

A still further aspect of the invention is an advanced video encoder comprising the method described above. This advanced video encoder may take the form of a specially designed video encoder signal processor embedded within a consumer device, such as a hand held high definition television (HDTV) device.

Yet another aspect of the invention is a method of intra mode prediction, that comprises:

(a) providing a 4×4 block of pixels [x] within a picture, wherein

${\lbrack x\rbrack \equiv \begin{bmatrix} X_{0,0} & X_{0,1} & X_{0,2} & X_{0,3} \\ X_{1,0} & X_{1,1} & X_{1,2} & X_{1,3} \\ X_{2,0} & X_{2,1} & X_{2,2} & X_{2,3} \\ X_{3,0} & X_{3,1} & X_{3,2} & X_{3,3} \end{bmatrix}};$

(b) providing 4 pixels

$\overset{\rightarrow}{h} = \begin{bmatrix} x_{0,{- 1}} \\ x_{1,{- 1}} \\ x_{2,{- 1}} \\ x_{3,{- 1}} \end{bmatrix}$

immediately to the left of the bock [x] in the picture when available;

(c) providing 4 pixels

$\overset{\rightarrow}{v} = \begin{bmatrix} x_{{- 1},0} \\ x_{{- 1},1} \\ x_{{- 1},2} \\ x_{{- 1},3} \end{bmatrix}$

immediately above the block [x] in the picture when available;

(d) providing a horizontal pixel position (H_(pos)) of the 4×4 block within the picture;

(e) providing a vertical pixel position (V_(pos)) of the 4×4 block within the picture; and

(f) outputting a lowest Sum of Absolute Transformed Differences (SATD) intra mode among intra modes while using [x], {right arrow over (h)}, {right arrow over (v)}, H_(pos), V_(pos), wherein the lowest SATD intra mode is selected from the group consisting of: i. a horizontal intra mode; ii. a vertical intra mode; and iii. a steady state (DC) intra mode.

The 4×4 pixel block [x] is typically oriented in the same sense as the original picture it is taken from, i.e. the indices increase in the same directions both horizontally and vertically. (Note that in the H.264 standard, the orientation is that pixel coordinate 0,0 is at the upper left of the picture.) In this orientation, the {right arrow over (h)} is immediately adjacent to the left side of the 4×4 pixel block [x], and the {right arrow over (v)} is immediately adjacent the top of the 4×4 pixel block [x].

It is readily apparent that the methods above may be easily modified for {right arrow over (h)} aligned on the right edge of the 4×4 pixel block [x], and the {right arrow over (v)} may be aligned with the bottom of the 4×4 pixel block [x]. Similarly, the picture or the 4×4 pixel block [x] indices may be readily reoriented with the orientation of pixel 0,0 being in any quadrant, while still performing the methods described here.

Additionally, the methods described herein may be readily applied to other H.264 pixel blocks, such as the 16×16 macroblock, or indeed generalized to block of size p×q where p and q may be independently selected from a group consisting of 2, 4, 8, 16, 32, 64, and 128.

In another aspect of the invention, the aforementioned steps may be incorporated into a more generalized method. The resultant generalized method of selecting an intra mode in a video image encoder may comprise:

(a) selecting a lowest Sum of Absolute Transformed Differences (SATD) intra mode among intra modes using inputs [x], H_(pos), V_(pos), {right arrow over (h)}, and {right arrow over (v)};

(b) wherein [x] is a p×q block of pixels within a picture with indices x_(i,j) for iε0, 1, 2, . . . , p−1, jε0, 1, 2, . . . , q−1; and

$\lbrack x\rbrack \equiv \begin{bmatrix} x_{0,0} & x_{0,1} & \ldots & x_{0,{q - 1}} \\ x_{1,0} & x_{1,1} & \ldots & x_{1,{q - 1}} \\ \vdots & \vdots & ⋰ & x_{2,{q - 1}} \\ x_{{p - 1},0} & x_{{p - 1},1} & x_{{p - 1},2} & x_{{p - 1},{q - 1}} \end{bmatrix}$

(c) wherein H_(pos) is a horizontal pixel position of the p×q block within the image;

(d) wherein V_(pos) is a vertical pixel position of the p×q block within the image;

(e) wherein {right arrow over (h)} is a horizontal vector immediately left of the p×q block [x], defined as {right arrow over (h)}≡(x_(0,−1), x_(1,−1), . . . , x_(p−1,−1)) relative to the indexing of the elements of [x];

(f) wherein {right arrow over (v)} is a horizontal vector immediately above the p×q block [x], defined as {right arrow over (v)}≡(x_(−1,0), x_(−1,1), . . . , x_(−1,q−1)) relative to the indexing of the elements of [x];

(g) wherein b is the bit depth of a pixel in the picture; (h) wherein the lowest SATD intra mode is determined among a group comprising:

-   -   i. a horizontal intra mode;     -   ii. a vertical intra mode; and     -   iii. a steady state (DC) intra mode; and

(i) outputting the lowest SATD intra mode to a computer readable medium.

The selecting step may comprise:

(a) calculating a horizontal predictor {right arrow over (H)}≡(H₀, H₁, . . . , H_(p−1))^(T), a vertical predictor {right arrow over (V)}≡(V₀, V₁, . . . , V_(q−1)), and a steady state (DC) predictor D;

(b) calculating a horizontal cost precursor C_(hs) and a vertical cost precursor C_(vs) using the horizontal predictor {right arrow over (H)}, the vertical predictor {right arrow over (V)}, and the steady state (DC) predictor D; and

(c) calculating a horizontal intra mode cost C_(H), a vertical intra mode cost C_(V), and a steady state (DC) intra mode cost C_(D) using the horizontal cost precursor C_(hs) and the vertical cost precursor C_(vs).

In the method above, the calculating a horizontal predictor {right arrow over (H)}, a vertical predictor {right arrow over (V)}, and a steady state (DC) predictor D may comprise:

(a) if H_(pos)≠0 and V_(pos)≠0 then: (i) setting {right arrow over (H)}=[T_(pq)]{right arrow over (h)} wherein {right arrow over (h)}≡(h₀, h₁, . . . , h_(p−1))^(T)≡(x_(0,0), x_(1,0), . . . , x_(p−1,0))^(T) (ii) setting {right arrow over (V)}=[T_(pq)]{right arrow over (v)} wherein {right arrow over (v)}≡(v₀, v₁, . . . , v_(q−1))≡(x_(0,0), x_(0,1), . . . , x_(0,q−1))^(T); and (iii) setting D=(H₀+V₀)/2;

(b) if H_(pos)=0 and V_(pos)≠0 then: (i) setting {right arrow over (H)}=(2¹⁵−1,0,0,0)^(T) wherein m is defined as 2^(m)≧p×p×q×2^(b+1); (ii) setting {right arrow over (V)}=[T_(pq)]{right arrow over (v)} wherein {right arrow over (v)}≡(v₀, v₁, . . . , v_(q))≡(x_(0,0), x_(0,1), . . . , x_(0,q))^(T); and (iii) setting D=V₀;

(c) if H_(pos)≠0 and V_(pos)=0 then: (i) setting {right arrow over (H)}=[T_(pq)]{right arrow over (h)} wherein {right arrow over (h)}≡(h₀, h₁, . . . , h_(p−1))^(T)≡(x_(0,0), x_(1,0), . . . , x_(p−1,0))^(T); (ii) setting {right arrow over (V)}=(2¹⁵−1,0,0,0)^(T) wherein n is defined as 2^(n)≧q×p×q×2^(b+1); and (iii) setting D=H₀; and

(d) if H_(pos)=0 and V_(pos)=0 then: (i) setting {right arrow over (H)}=(2¹⁵−1,0,0,0)^(T) wherein m is defined as 2^(m)≧p×p×q×2^(b+1); (ii) setting {right arrow over (V)}=(2¹⁵−1,0,0,0)^(T) wherein n is defined as 2^(n)≧q×p×q×2^(b+1); and (iii) setting D=128×16.

The horizontal cost precursor C_(hs) and the vertical cost precursor C_(vs) may be calculated by steps comprising:

-   -   (a) calculating the values X_(i,0), X_(0,j) for iε0, 1, 2, . . .         , p−1, jε0, 1, 2, . . . , q−1 using the relationship         [X]=[T_(pq)][x][T_(pq)]^(T) where

${\lbrack X\rbrack \equiv \begin{bmatrix} X_{0,0} & X_{0,1} & \ldots & X_{0,{q - 1}} \\ X_{1,0} & X_{1,1} & \ldots & X_{1,{q - 1}} \\ \vdots & \vdots & ⋰ & X_{2,{q - 1}} \\ X_{{p - 1},0} & X_{{p - 1},1} & X_{{p - 1},2} & X_{{p - 1},{q - 1}} \end{bmatrix}};$

-   -   (b) calculating the horizontal cost precursor

${C_{hs} = {\sum\limits_{i = 1}^{p - 1}{X_{i,0}}}};$

and

-   -   (c) calculating the vertical cost precursor

$C_{vs} = {\sum\limits_{j = 1}^{q - 1}{{X_{0,j}}.}}$

The horizontal intra mode cost C_(H) may be calculated using steps comprising: calculating

$C_{H} = {{\sum\limits_{i = 0}^{p - 1}{{H_{i} - X_{i,0}}}} + {C_{vs}.}}$

Similarly, the vertical intra mode cost C_(V) may be calculated with steps comprising: calculating

$C_{v} = {{\sum\limits_{j = 0}^{q - 1}{{V_{j} - X_{0,j}}}} + {C_{hs}.}}$

Finally, the steady state (DC) intra mode cost C_(D) may be calculated with steps comprising: calculating C_(D)=|D−X_(0,0)|+C_(hs)+C_(vs).

In the method above, the calculating means may comprise: selecting the lowest SATD intra mode with a lowest associated intra mode cost among the group consisting of: the horizontal intra mode cost C_(H), the vertical intra mode cost C_(V), and the steady state (DC) intra mode cost C_(D).

In another aspect of the invention, a computer readable medium comprising a programming executable may be capable of performing on a computer the methods described above.

In still another aspect of the invention, the methods of intra mode prediction above may be generalized for: (a) where the p dimension is selected from a group of dimensions consisting of: 2, 4, 8, 16, 32, 64, and 128; and (b) where the q dimension is selected from a group of dimensions consisting of: 2, 4, 8, 16, 32, 64, and 128. Still higher dimensions are also possible, however, they are likely computationally prohibitive nature of the increasingly complex calculations.

Further aspects of the invention will be brought out in the following portions of the specification, wherein the detailed description is for the purpose of fully disclosing preferred embodiments of the invention without placing limitations thereon.

BRIEF DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWING(S)

The invention will be more fully understood by reference to the following drawings which are for illustrative purposes only:

FIG. 1A is a flow chart of a Hadamard transform of a set of 4 pixels into a Hadamard transform of the 4 pixels.

FIG. 1B is a flow chart of a pruned Hadamard transform of a 4×4 block of pixels that produces only the vertical column and horizontal row components of the Hadamard transform of the 4×4 block of pixels.

FIG. 2A is an overview flow chart of the intra mode selection method.

FIG. 2B is a flow chart of the selected intra mode being output to one or more computer readable media.

FIG. 3A is a matrix of the 4×1 vector {right arrow over (h)} to the left to the 4×4 block and 1×4 element vector {right arrow over (v)} above the 4×4 block.

FIG. 3B is a matrix of the left transform coefficients and top transform coefficients that correspond to the left 4×1 and top 1×4 elements of FIG. 3A depicting the relationship between the spatial and frequency domain intra predictor for the horizontal and vertical modes.

FIG. 4 is a flowchart that details the computation of the frequency domain predictors for the intra vertical, horizontal, and steady state (or DC) intra modes.

FIG. 5 is a flowchart that predicts the SATD costs of the various horizontal, vertical, or DC predictions 500. Using these costs, it outputs a selected intra mode with the least SATD.

DETAILED DESCRIPTION OF THE INVENTION

Referring more specifically to the drawings, for illustrative purposes the present invention is embodied in the method generally depicted in FIG. 1A through FIG. 5. It will be appreciated that the apparatus may vary as to configuration and as to details of the parts, and that the method may vary as to the specific steps and sequence, without departing from the basic concepts as disclosed herein.

Introduction

In the H.264/AVC video compression standard, the 4×4 intra mode decision is made by performing Lagrangian rate-distortion optimization on each of the 9 intra mode possibilities. The intra mode selected is the one that provides the lowest distortion (SATD) between the original 4×4 input block, and the reconstructed block for a particular intra mode. This process is extremely computationally expensive.

The process described below reduces the computational cost of the 4×4 intra mode decision significantly, from 24 1D 4 point Hadamard transforms to merely 4 of them.

DEFINITIONS

“Computer” means any device capable of performing the steps, methods, or producing signals as described herein, including but not limited to: a microprocessor, a microcontroller, a video processor, a digital state machine, a field programmable gate array (FGPA), a digital signal processor, a collocated integrated memory system with microprocessor and analog or digital output device, a distributed memory system with microprocessor and analog or digital output device connected by digital or analog signal protocols.

“Computer readable medium” means any source of organized information that may be processed by a computer to perform the steps described herein to result in, store, perform logical operations upon, or transmit, a flow or a signal flow, including but not limited to: random access memory (RAM), read only memory (ROM), a magnetically readable storage system; optically readable storage media such as punch cards or printed matter readable by direct methods or methods of optical character recognition; other optical storage media such as a compact disc (CD), a digital versatile disc (DVD), a rewritable CD and/or DVD; electrically readable media such as programmable read only memories (PROMs), electrically erasable programmable read only memories (EEPROMs), field programmable gate arrays (FGPAs), flash random access memory (flash RAM); and information transmitted by electromagnetic or optical methods including, but not limited to, wireless transmission, copper wires, and optical fibers.

“SATD” means the Sum of Absolute Transformed Differences, which is a widely used video quality metric used for block-matching in motion estimation for video compression. It works by taking a frequency transform, usually a Hadamard transform, of the differences between the pixels in the original block and the corresponding pixels in the block being used for comparison. The transform itself is often of a small block rather than the entire macroblock. For example, in H.264/AVC, a series of 4×4 blocks are transformed rather than doing a more processor-intensive 16×16 transform.

“DCT” (Discrete Cosine Transformation) means a process that converts images from three-dimensions (3D) to two-dimensions (2D) by using the Discrete Cosine (DC) coefficient to examine the luminance of each block of pixels used to form an image. This process is typically used in MPEG and JPEG image compression.

“Hadamard Transforms” are defined below. Although the Hadamard transform has numerous applications to signal processing and analysis, this application describes specific applications of the Hadamard Transform in the context of video processing, without limitation of other, more general, applications.

Hadamard transforms are used in the intra mode estimation in the Joint Model (JM) of the Advance Video Coding (AVC) encoder. Refer now to FIG. 1A, a Hadamard transform takes an input of 4 pixels (otherwise referred to as pels) and transforms them into a Hadamard Transform of 4 pixels. Similarly, in FIG. 1B, a 4×4 block may also be selectively transformed into a Hadamard Transform of a 4×4 block while outputting only a particular horizontal row and vertical column. In this instance, selected rows or columns of the input 4×4 block may be transformed as described further below. The Hadamard transform inherently possesses characteristics of uniform scaling such that each transform coefficient has the same amplification.

The Hadamard Vector Transform of a 4 pixel input may mathematically be defined in the following manner. Let {right arrow over (s)}=[s₀,s₁,s₂,s₃]^(T) be a vector consisting of 4 elements. The Hadamard transform [T₄] of {right arrow over (s)} is defined as {right arrow over (S)}=[S₀,S₁,S₂,S₃]^(T)=[T₄]{right arrow over (s)} where

$\left\lbrack T_{4} \right\rbrack = {\begin{bmatrix} 1 & 1 & 1 & 1 \\ 1 & 1 & {- 1} & {- 1} \\ 1 & {- 1} & {- 1} & 1 \\ 1 & {- 1} & 1 & {- 1} \end{bmatrix}.}$

The Hadamard Transformation of a 4×4 Block is defined as follows.

Let [y] be a 4×4 block of pixels such that

$\lbrack y\rbrack = {\begin{bmatrix} y_{0,0} & y_{0,1} & y_{0,2} & y_{0,3} \\ y_{1,0} & \; & \; & \; \\ y_{2,0} & \; & ⋰ & \; \\ y_{3,0} & \; & \; & y_{3,3} \end{bmatrix}.}$

The Hadamard transform H_(4×4)(y) of [y] is defined as

$\lbrack Y\rbrack = {\begin{bmatrix} Y_{0,0} & Y_{0,1} & Y_{0,2} & Y_{0,3} \\ Y_{1,0} & \; & \; & \; \\ Y_{2,0} & \; & ⋰ & \; \\ Y_{3,0} & \; & \; & Y_{3,3} \end{bmatrix} = {{T_{4 \times 4}\left( \lbrack y\rbrack \right)} = {{\left\lbrack T_{4} \right\rbrack \lbrack y\rbrack}\left\lbrack T_{4} \right\rbrack}^{T}}}$

where

${\begin{bmatrix} Y_{0,0} \\ Y_{1,0} \\ Y_{2,0} \\ Y_{3,0} \end{bmatrix} = {{\left\lbrack T_{4} \right\rbrack \lbrack y\rbrack}\overset{\rightarrow}{u}}},{\begin{bmatrix} Y_{0,0} & Y_{0,1} & Y_{0,2} & Y_{0,3} \end{bmatrix} = {{{\overset{\rightarrow}{u}}^{T}\lbrack y\rbrack}\left\lbrack T_{4} \right\rbrack}^{T}},$

and {right arrow over (u)}^(T)=[1 1 1 1].

Referring back to FIG. 1B, to selectively determine row zero and column zero of the Hadamard block transform. Here, an input 4×4 pixel block, such as [y], is input. A selection vector {right arrow over (u)} is used to selectively determine the column zero Hadamard transform of the input 4×4 pixel block [y], to produce

$\begin{bmatrix} Y_{0,0} \\ Y_{1,0} \\ Y_{2,0} \\ Y_{3,0} \end{bmatrix}.$

Similarly, selection vector {right arrow over (u)}^(T) is used with the input 4×4 pixel block [y], to produce the row zero Hadamard transform of the input 4×4 pixel block [y], which results in [Y_(0,0) Y_(0,1) Y_(0,2) Y_(0,3)].

The combination of the row zero block Hadamard transform, and the column zero block Hadamard transform is denoted as T_(RC), indicating that only a row and column of the block transformation is to be processed. This may also be referred to as a “pruned” Hadamard transformation.

Process Overview

Referring now to FIG. 2, an overview of an embodiment of a method of determining an optimal intra mode selection 200 is shown. In the embodiment shown, an input 4×4 block of pixels 202 is used as an input to the pruned 4×4 Hadamard block transformation 204 to produce the vertical and horizontal frequency component of the 4×4 Hadamard block transform output 206. This output 206 will be used subsequently as described below.

The 4×1 pixels 208 (the 4 top elements immediately above the input 4×4 block of pixels 202) are used as input into a T₄ Hadamard transform 210 to produce a vertical prediction Hadamard transform output 212.

Similarly, the left 1×4 pixels 214 (the 4 left elements immediately left of the input 4×4 block of pixels 202) are used as input into a T₄ Hadamard transform 216 to produce a horizontal prediction Hadamard transform output 218.

T₄ Hadamard vertical 212 and horizontal 218 predictions are used to estimate the DC prediction 222.

The following inputs are compared 224 to determine the optimal intra mode prediction 226: 1) the pruned 4×4 Hadamard block transform output 206; 2) the vertical prediction Hadamard transform output 212; 3) the horizontal prediction Hadamard transform output 218; and 4) the DC prediction 222.

Only horizontal 218, vertical 212, and DC 222 predictions are used in the intra DCT mode decision coefficients. The intra predictions are performed in the frequency domain.

Refer now to FIG. 2B, which is a diagram of the various computer readable media that the selected intra mode 226 may be output. It may be sent to a microwave transmitter 228 for propagation through microwaves 230 to a distant connection. The intra mode 226 may be stored on a memory device faster or slower than the depicted USB storage device 232, as well as stored on a digital or analog recording medium exemplified as a high definition video recorder 234, or a still or moving digital camera 236. Similarly, the intra mode 226 may be displayed, with or without further processing, on a graphical device such as a plasma or liquid crystal display 238. The intra mode 226 may be transmitted through a network 240, which may be wireless, or may use TCP/IP or other interconnection techniques to any other device in this FIG. 2B. Still another output of the intra mode 226 may be to a multiple media device 242, which may have a VHS video recorder 244, or a DVD recorder 246. Yet another output of the intra mode 226 may be to a computer 248 that may or may not be connected to a display device 250, to be ultimately viewed by a person 252. The computer 248 may itself have a wide variety of computer readable media available for subsequent secondary output of the intra mode 226. Finally, the intra mode 226 may be transmitted through a wired connection 254, which may have one or more conductors 256 of electromagnetic information, ultimately being transmitted to one or more of the devices present in this FIG. 2B.

Estimate the Intra Macroblock DCT Coefficients

To reduce computation, only horizontal, vertical, and DC predictions are used in the estimation of the intra DCT coefficients. In particular, the intra predictions are computed in frequency domain; the DC prediction is derived from the horizontal and vertical predictions. And, finally, the prediction residue with the minimal SATD is selected as the intra mode.

Refer now to FIGS. 3A and 3B, which taken together describe the relationship 300 between the spatial and frequency domain intra predictor for horizontal and vertical modes. Here, an initial spatial domain representation (in FIG. 3A) of an input 4×4 block 302 is shown as [x] with elements x_(i,j), where i,jε(0,1,2,3). The frequency domain representation (in FIG. 3B) of the 4×4 transformation 304 is shown as the transformed matrix[X], with elements X_(i,j), where i,jε(0,1,2,3).

For convenience, the 4×1 column vector to the left of vector (x_(0,0),x_(1,0),x_(2,0),x_(3,0))^(T) with elements (x_(0,−1),x_(1,−1),x_(2,−1),x_(3,−1))^(T) is denoted as {right arrow over (h)}=(h₀,h₁,h₂,h₃)^(T) 306. The Hadamard transform of h is denoted as {right arrow over (H)}=(H₀,H₁,H₂,H₃)^(T) 308.

Similarly, the 1×4 row vector above 4×4 block [x] 302 are (x_(−1,0),x_(−1,1),x_(−1,2),x_(−1,3)), which are for convenience denoted 310 as {right arrow over (v)}=(v₀,v₁,v₂,v₃). The Hadamard transform of {right arrow over (v)} 312 is denoted {right arrow over (V)}=(V₀,V₁,V₂,V₃)^(T).

Compute Frequency Domain Predictors for the Intra Vertical, Horizontal, and DC Prediction Modes

This process may be followed more readily by referring to FIG. 4, which details the computation of the frequency domain predictors for the intra vertical, horizontal, and steady state (or DC) modes 400.

First, input scalar index positions (H_(pos),V_(pos)) of the top left pixels of a 4×4 pixel block [x] in a picture that begins with pixels 0,0 (the upper left corner of the picture in the H.264 design specification) and continues to pixel position values m,n. Also input the pixel block [x] 402.

Next, from the 4 pixels immediately to the left and above the 4×4 pixel block [x] are denoted as 404 {right arrow over (h)}=(h₀,h₁, h₂, h₃)^(T), {right arrow over (v)}=(v₀,v₁,v₂,v₃).

At this point, now calculate the horizontal predictor {right arrow over (H)}=[H₀,H₁,H₂,H₃]^(T), the vertical predictor {right arrow over (V)}=[V₀,V₁,V₂,V₃]^(T), and the steady state (DC) predictor D as follows:

If H_(pos)≠0 (406) and V_(pos)≠0 (408), then:

{right arrow over (H)}=[T₄]{right arrow over (h)}

{right arrow over (V)}=[T₄]{right arrow over (v)}

D=(H ₀ +V ₀)/2

at (410).

If H_(pos)=0 (e.g., not H_(pos)≠0 at 406) and V_(pos)≠0 (at 412), then:

{right arrow over (H)}=[2¹⁵−1,0,0,0]^(T)

{right arrow over (V)}=[T₄]{right arrow over (v)}

D=V₀

at (414).

If H_(pos)≠0 (406) and V_(pos)=0 (e.g., not V_(pos)≠0 at 408), then:

{right arrow over (H)}=[T₄]{right arrow over (h)}

{right arrow over (V)}=[2¹⁵−1,0,0,0]^(T)

D=H₀

at (416).

If H_(pos)=0 (e.g., not H_(pos)≠0 at 406) and V_(pos)=0 (e.g., not V_(pos)≠0 at 412), then:

{right arrow over (H)}=[2¹⁵−1,0,0,0]^(T)

{right arrow over (V)}=[2¹⁵−1,0,0,0]^(T)

D=128×16

at (418).

Here, it is assumed that the pixels can only take on 8 bits of information. In particular, the DC predictor 2¹⁵−1 appearing in block 418 corresponds to the DC prediction for 8 bits per pixel. The predictor {right arrow over (H)}=[2¹⁵−1,0,0,0]^(T) in 414, 418, and the predictor {right arrow over (V)}=[2¹⁵−1,0,0,0]^(T) in 416, 418, are selected to make sure that they will have sufficiently large intra prediction cost for 8 bits per pixel, and consequently the corresponding prediction mode will not be selected as the minimal cost intra prediction mode in FIG. 5. This is consistent with the H.264/AVC standard.

Regardless of which calculation branch was taken from 410, 414, 416, or 418, next the cost is calculated 420.

Compute Intra Prediction Cost

Refer now to FIG. 5, which predicts the computational costs of the various horizontal, vertical, or DC predictions 500, and using these, outputs a selected intra mode with least SATD. To this evaluation is provided the {right arrow over (H)}, {right arrow over (V)}, D values determined above, as well as the input 4×4 pixel block [x]. Next, the values of X_(i,0), X_(0,i) for iε0, 1, 2, 3 are determined 504 using the relationship described above for the pruned Hadamard block transform, i.e.

$\begin{bmatrix} X_{0,0} \\ X_{1,0} \\ X_{2,0} \\ X_{3,0} \end{bmatrix} = {{{\left\lbrack T_{4} \right\rbrack \lbrack x\rbrack}\overset{\rightarrow}{u}\mspace{14mu} {{and}\mspace{14mu}\begin{bmatrix} X_{0,0} & X_{0,1} & X_{0,2} & X_{0,3} \end{bmatrix}}} = {{{\overset{\rightarrow}{u}\lbrack x\rbrack}\left\lbrack T_{4} \right\rbrack}^{T}.}}$

Cost precursors are then formed 506

$C_{hs} = {{\sum\limits_{i = 1}^{3}{{X_{i,0}}\mspace{14mu} {and}\mspace{14mu} C_{vs}}} = {\sum\limits_{j = 1}^{3}{{X_{0,j}}.}}}$

Finally, the costs are calculated 508, where the cost of the horizontal prediction is

${C_{H} = {{\sum\limits_{i = 0}^{3}\; {{H_{i} - X_{i,0}}}} + C_{vs}}},$

the cost of the vertical prediction is

${C_{v} = {{\sum\limits_{j = 0}^{3}\; {{V_{j} - X_{0,j}}}} + C_{hs}}},$

and the cost of the DC prediction is C_(D)=[D−X_(0,0)]+C_(hs)+C_(vs).

Once the predicted costs are determine, the appropriate intra mode is selected from the group of Horizontal Prediction, Vertical Prediction, and DC Prediction.

Minimal Cost Intra Mode Prediction

The intra mode with the minimal cost is selected as the intra prediction mode. In particular:

If C_(H)≦C_(V) and C_(H)≦C_(D), then select Horizontal Prediction 510;

If C_(H)≦C_(V) and C_(H)>C_(D), then select DC Prediction 512;

If C_(H)>C_(V) and C_(V)≦C_(D), then select Vertical Prediction 514; and finally,

If C_(H)>C_(V) and C_(V)>C_(D), select DC Prediction 516.

The minimal cost prediction among C_(H), C_(V), and C_(D) is then output as the appropriate associated predicted intra mode. From this point, the selected intra mode is used within the advanced video encoder to compress the 4×4 block.

Discussion

As previously described, the H.264/AVC video compression standard the 4×4 intra mode decision is made by performing Lagrangian rate-distortion optimization on each of the 9 intra mode possibilities. The intra mode selected is the one that provides the lowest SATD between the original 4×4 input block, and the reconstructed block for a particular intra mode, where SATD is directly computed as the sum of absolute value of the transform of the difference. This process is extremely computationally expensive.

For selecting among the horizontal, vertical, and DC intra modes with the direct method, each of SATD of these intra modes is separately computed with a 4×4 Hadamard block transform for a total of three 2D Hadamard transforms.

Each 2D 4×4 Hadamard transform is computationally equivalent to 8 of the 1D 4 point Hadamard transforms when row and column decomposition is used for the 2d transform.

In this process, the SATD is computed as the sum of the absolute value of difference of the transform coefficients, and only the horizontal, vertical, and DC transform coefficients are computed using 1D 4 point Hadamard transform. Thus, this process reduces the number of equivalent 1D Hadamard transforms of the horizontal, vertical, and DC intra modes from 24 to 4 by performing the intra mode prediction in the frequency domain, while providing the same intra mode selection as direct calculations of SATD.

CONCLUSION

Although the description above contains many details, these should not be construed as limiting the scope of the invention but as merely providing illustrations of some of the presently preferred embodiments of this invention. Therefore, it will be appreciated that the scope of the present invention fully encompasses other embodiments which may become obvious to those skilled in the art. In the appended claims, reference to an element in the singular is not intended to mean “one and only one” unless explicitly so stated, but rather “one or more.” All equivalents to the elements of the above-described preferred embodiments that are known to those of ordinary skill in the art are expressly incorporated herein by reference and are intended to be encompassed by the present claims. Moreover, it is not necessary for a device or method to address each and every problem sought to be solved by the present invention, for it to be encompassed by the present claims. Furthermore, no element, component, or method step in the present disclosure is intended to be dedicated to the public regardless of whether the element, component, or method step is explicitly recited in the claims. No claim element herein is to be construed under the provisions of 35 U.S.C. 112, sixth paragraph, unless the element is expressly recited using the phrase “means for.” 

1. A method of selecting an intra mode in a video image encoder, comprising: (a) selecting a lowest Sum of Absolute Transformed Differences (SATD) intra mode among intra modes using inputs [x], H_(pos), V_(pos), {right arrow over (h)}, and {right arrow over (v)}; (b) wherein [x] is a 4×4 block of pixels within a picture and $\lbrack x\rbrack \equiv \begin{bmatrix} x_{0,0} & x_{0,1} & x_{0,2} & x_{0,3} \\ x_{1,0} & x_{1,1} & x_{1,2} & x_{1,3} \\ x_{2,0} & x_{2,1} & x_{2,2} & x_{2,3} \\ x_{3,0} & x_{3,1} & x_{3,2} & x_{3,3} \end{bmatrix}$ (c) wherein H_(pos) is a horizontal pixel position of the 4×4 block within the image; (d) wherein V_(pos) is a vertical pixel position of the 4×4 block within the image; (e) wherein {right arrow over (h)} is a horizontal vector immediately left of the 4×4 block [x], defined as {right arrow over (h)}≡(x_(0,−1),x_(1,−1),x_(2,−1),x_(3,−1))^(T) relative to the indexing of the elements of [x]; (f) wherein {right arrow over (h)} is a horizontal vector immediately above the 4×4 block [x], defined as {right arrow over (v)}≡(x_(−1,0),x_(−1,1),x_(−1,2),x_(−1,3))^(T) relative to the indexing of the elements of [x]; and (g) wherein the lowest SATD intra mode is determined among a group comprising: i. a horizontal intra mode; ii. a vertical intra mode; and iii. a steady state (DC) intra mode; and (h) outputting the lowest SATD intra mode to a computer readable medium.
 2. The method of intra mode prediction of claim 1, wherein said outputting step comprises: (a) calculating a horizontal predictor {right arrow over (H)}≡(H₀,H₁,H₂,H₃)^(T), a vertical predictor {right arrow over (V)}≡(V₀,V₁,V₂,V₃)^(T), and a steady state (DC) predictor D; (b) calculating a horizontal cost precursor C_(hs) and a vertical cost precursor C_(vs) using the horizontal predictor {right arrow over (H)}, the vertical predictor {right arrow over (V)}, and the steady state (DC) predictor D; and (c) calculating a horizontal intra mode cost C_(H), a vertical intra mode cost C_(V), and a steady state (DC) intra mode cost C_(D) using the horizontal cost precursor C_(hs) and the vertical cost precursor C_(vs).
 3. The method of intra mode prediction of claim 2, wherein the calculating of the horizontal predictor {right arrow over (H)}, the vertical predictor {right arrow over (V)}, and the steady state (DC) predictor D comprises: if H_(pos)≠0 and V_(pos)≠0 then: i. setting {right arrow over (H)}≡(H₀,H₁,H₂,H₃)^(T)=[T₄]{right arrow over (h)} where {right arrow over (h)}≡(h₀,h₁,h₂,h₃)^(T)≡(x_(0,−1),x_(1,−1),x_(2,−1),x_(3,−1))^(T); ii. setting {right arrow over (V)}≡(V₀,V₁,V₂,V₃)^(T)=[T₄]{right arrow over (v)} where {right arrow over (v)}≡(v₀,v₁,v₂,v₃)^(T)≡(x_(−1,0),x_(−1,1),x_(−1,2),x_(−1,3))^(T); and iii. setting D=(H₀+V₀)/2.
 4. The method of intra mode prediction of claim 2, wherein the calculating of the horizontal predictor {right arrow over (H)}, the vertical predictor {right arrow over (V)}, and the steady state (DC) predictor D comprises: if H_(pos)=0 and V_(pos)≠0 then: i. setting {right arrow over (H)}=(2¹⁵−1,0,0,0)^(T); ii. setting {right arrow over (V)}≡(V₀,V₁,V₂,V₃)^(T)=[T₄]{right arrow over (v)} where {right arrow over (v)}≡(v₀,v₁,v₂,v₃)^(T)≡(x_(−1,0),x_(−1,1),x_(−1,2),x_(−1,3))^(T); and iii. setting D=V₀.
 5. The method of intra mode prediction of claim 2, wherein the calculating of the horizontal predictor {right arrow over (H)}, the vertical predictor {right arrow over (V)}, and the steady state (DC) predictor D comprises: if H_(pos)≠0 and V_(pos)=0 then: i. setting {right arrow over (H)}≡(H₀,H₁,H₂,H₃)^(T)=[T₄]{right arrow over (h)} where {right arrow over (h)}≡(h₀,h₁,h₂,h₃)^(T)≡(x_(0,−1),x_(1,−1),x_(2,−1),x_(3,−1))^(T); ii. setting {right arrow over (V)}=(2¹⁵−1,0,0,0)^(T); and iii. setting D=H₀.
 6. The method of intra mode prediction of claim 2, wherein the calculating of the horizontal predictor {right arrow over (H)}, the vertical predictor {right arrow over (V)}, and v steady state (DC) predictor D comprises: if H_(pos)=0 and V_(pos)=0 then: i. setting {right arrow over (H)}=(2¹⁵−1,0,0,0)^(T); ii. setting {right arrow over (V)}=(2¹⁵−1,0,0,0)^(T); and iii. setting D=128×16.
 7. The method of intra mode prediction of claim 2, wherein the calculating of the horizontal cost precursor C_(hs) and the vertical cost precursor C_(vs) comprises: (a) calculating the values X_(i,0), X_(0,i) for iε0, 1, 2, 3 using the relationships ${\begin{bmatrix} X_{0,0} \\ X_{1,0} \\ X_{2,0} \\ X_{3,0} \end{bmatrix} = {{{\left\lbrack T_{4} \right\rbrack \lbrack x\rbrack}\overset{\rightarrow}{u}\mspace{14mu} {{and}\mspace{14mu}\begin{bmatrix} X_{0,0} & X_{0,1} & X_{0,2} & X_{0,3} \end{bmatrix}}} = {{\overset{\rightarrow}{u}\lbrack x\rbrack}\left\lbrack T_{4} \right\rbrack}^{T}}},$ (i) wherein {right arrow over (u)}^(T)=[1 1 1 1], (ii) wherein ${\left\lbrack T_{4} \right\rbrack = \begin{bmatrix} 1 & 1 & 1 & 1 \\ 1 & 1 & {- 1} & {- 1} \\ 1 & {- 1} & {- 1} & 1 \\ 1 & {- 1} & 1 & {- 1} \end{bmatrix}};$ (b) calculating the horizontal cost precursor ${C_{hs} = {\sum\limits_{i = 1}^{3}\; {X_{i,0}}}};$  and (c) calculating the vertical cost precursor $C_{vs} = {\sum\limits_{j = 1}^{3}\; {{X_{0,j}}.}}$
 8. The method of intra mode prediction of claim 2, wherein the calculating the horizontal intra mode cost C_(H) comprises: calculating $C_{H} = {{\sum\limits_{i = 0}^{3}{{H_{i} - X_{i,0}}}} + {C_{vs}.}}$
 9. The method of intra mode prediction of claim 2, wherein the calculating the vertical intra mode cost C_(V) comprises: calculating $C_{v} = {{\sum\limits_{j = 0}^{3}\; {{V_{j} - X_{0,j}}}} + {C_{hs}.}}$
 10. The method of intra mode prediction of claim 2, wherein the calculating the steady state (DC) intra mode cost C_(D) comprises: calculating C_(D)=|D−X_(0,0)|+C_(hs)+C_(vs).
 11. The method of intra mode prediction of claim 2, further comprising: selecting the lowest SATD intra mode with a lowest associated intra mode cost among the group consisting of: the horizontal intra mode cost C_(H), the vertical intra mode cost C_(V), and the steady state (DC) intra mode cost C_(D).
 12. A computer readable medium comprising programming executable on a computer and capable of causing said computer to perform the method of claim
 1. 13. An advanced video encoder comprising a computer and programming executable by said computer for carrying out the method of claim
 1. 14. A method of selecting an intra mode in a video image encoder, comprising: (a) selecting a lowest Sum of Absolute Transformed Differences (SATD) intra mode among intra modes using inputs [x], H_(pos), V_(pos), {right arrow over (h)}, and {right arrow over (v)}; (b) wherein [x] is a p×q block of pixels within a picture with indices x_(i,j) for iε0, 1, 2, . . . , p−1, jε0, 1, 2, . . . , q−1; and $\lbrack x\rbrack \equiv \begin{bmatrix} x_{0,0} & x_{0,1} & \cdots & x_{0,{q - 1}} \\ x_{1,0} & x_{1,1} & \cdots & x_{1,{q - 1}} \\ \vdots & \vdots & ⋰ & x_{2,{q - 1}} \\ x_{{p - 1},0} & x_{{p - 1},1} & x_{{p - 1},2} & x_{{p - 1},{q - 1}} \end{bmatrix}$ (c) wherein H_(pos) is a horizontal pixel position of the p×q block within the image; (d) wherein V_(pos) is a vertical pixel position of the p×q block within the image; (e) wherein {right arrow over (h)} is a horizontal vector immediately left of the p×q block [x], defined as {right arrow over (h)}≡(x_(0,−1), x_(1,−1), . . . , x_(p−1,−1))^(T) relative to the indexing of the elements of [x]; (f) wherein {right arrow over (v)} is a horizontal vector immediately above the p×q block [x], defined as {right arrow over (v)}≡(x_(−1,0), x_(−1,1), . . . , x_(−1,q−1))^(T) relative to the indexing of the elements of [x]; (g) wherein b is the bit depth of a pixel in the picture; (h) wherein the lowest SATD intra mode is determined among a group comprising: i. a horizontal intra mode; ii. a vertical intra mode; and iii. a steady state (DC) intra mode; and (i) outputting the lowest SATD intra mode to a computer readable medium.
 15. The method of intra mode prediction of claim 14, wherein the selecting step comprises: (a) calculating a horizontal predictor {right arrow over (H)}≡(H₀, H₁, . . . , H_(p−1))^(T), a vertical predictor {right arrow over (V)}≡(V₀, V₁, . . . , V_(q−1))^(T), and a steady state (DC) predictor D; (b) calculating a horizontal cost precursor C_(hs) and a vertical cost precursor C_(vs) using the horizontal predictor {right arrow over (H)}, the vertical predictor {right arrow over (V)}, and the steady state (DC) predictor D; and (c) calculating a horizontal intra mode cost C_(H), a vertical intra mode cost C_(V), and a steady state (DC) intra mode cost C_(D) using the horizontal cost precursor C_(hs) and the vertical cost precursor C_(vs).
 16. The method of intra mode prediction of claim 15, wherein the calculating a horizontal predictor {right arrow over (H)}, a vertical predictor {right arrow over (V)}, and a steady state (DC) predictor D comprise: (a) if H_(pos)≠0 and V_(pos)≠0 then: i. setting {right arrow over (H)}=[T_(pq)]{right arrow over (h)} wherein {right arrow over (h)}≡(h₀, h₁, . . . , h_(p−1))^(T)≡(x_(0,0), x_(1,0), . . . , x_(p−1,0))^(T) ii. setting {right arrow over (V)}=[T_(pq)]{right arrow over (v)} wherein {right arrow over (v)}≡(v₀, v₁, . . . , v_(q−1))≡(x_(0,0), x_(0,1), . . . , x_(0,q−1))^(T); and iii. setting D=(H₀+V₀)/2; (b) if H_(pos)=00 and V_(pos)≠0 then: i. setting {right arrow over (H)}=(2¹⁵−1,0,0,0)^(T) wherein m is defined as 2^(m)≧p×p×q×2^(b+1); ii. setting {right arrow over (V)}=[T_(pq)]{right arrow over (v)} wherein {right arrow over (v)}≡(v₀, v₁, . . . , v_(q))≡(x_(0,0), x_(0,1), . . . , x_(0,q))^(T); and iii. setting D=V₀; (c) if H_(pos)≠0 and V_(pos)=0 then: i. setting {right arrow over (H)}=[T_(pq)]{right arrow over (h)} wherein {right arrow over (h)}≡(h₀, h₁, . . . , h_(p−1))^(T)≡(x_(0,0), x_(1,0), . . . , x_(p−1,0))^(T); ii. setting {right arrow over (V)}=(2¹⁵−1,0,0,0)^(T) wherein n is defined as 2^(n)>q×p×q×2^(b+1); and iii. setting D=H₀; (d) if H_(pos)=0 and V_(pos)=0 then: i. setting {right arrow over (H)}=(2¹⁵−1,0,0,0)^(T) wherein m is defined as 2^(m)>p×p×q'2^(b+1); ii setting {right arrow over (V)}=(2¹⁵−1,0,0,0)^(T) wherein n is defined as 2^(n)>q×p×q×2^(b+1); and iii. setting D=128×16.
 17. The method of intra mode prediction of claim 15, wherein the calculating the horizontal cost precursor C_(hs) and the vertical cost precursor C_(vs) comprises: (a) calculating the values X_(i,0), X_(0,j) for iε(0, 1, 2, . . . , p−1, jε0, 1, 2, . . . , q−1 using the relationship [X]=[T_(pq)][x][T_(pq)]^(T), where ${\lbrack X\rbrack \equiv \begin{bmatrix} X_{0,0} & X_{0,1} & \cdots & X_{0,{q - 1}} \\ X_{1,0} & X_{1,1} & \cdots & X_{1,{q - 1}} \\ \vdots & \vdots & ⋰ & X_{2,{q - 1}} \\ X_{{p - 1},0} & X_{{p - 1},1} & X_{{p - 1},2} & X_{{p - 1},{q - 1}} \end{bmatrix}};$ (b) calculating the horizontal cost precursor ${C_{hs} = {\sum\limits_{i = 1}^{p - 1}\; {X_{i,0}}}};$  and (c) calculating the vertical cost precursor $C_{vs} = {\sum\limits_{j = 1}^{q - 1}\; {{X_{0,j}}.}}$
 18. The method of intra mode prediction of claim 15, wherein the calculating the horizontal intra mode cost C_(H) comprises: calculating $C_{H} = {{\sum\limits_{i = 0}^{p - 1}\; {{H_{i} - X_{i,0}}}} + {C_{vs}.}}$
 19. The method of intra mode prediction of claim 15, wherein the calculating the vertical intra mode cost C_(V) comprises: calculating $C_{v} = {{\sum\limits_{j = 0}^{q - 1}\; {{V_{j} - X_{0,j}}}} + {C_{hs}.}}$
 20. The method of intra mode prediction of claim 15, wherein the calculating the steady state (DC) intra mode cost C_(D) comprises: calculating C_(D)=D−X_(0,0)+C_(hs)+C_(vs).
 21. The method of intra mode prediction of claim 15, wherein the calculating means comprises: selecting the lowest SATD intra mode with a lowest associated intra mode cost among the group consisting of: the horizontal intra mode cost C_(H), the vertical intra mode cost C_(V), and the steady state (DC) intra mode cost C_(D).
 22. A computer readable medium comprising a programming executable capable of performing on a computer the method of claim
 14. 23. The method of intra mode prediction of claim 15, wherein: (a) the p dimension is selected from a group of dimensions consisting of: 2, 4, 8, 16, 32, 64, and 128; and (a) the q dimension is selected from a group of dimensions consisting of: 2, 4, 8, 16, 32, 64, and
 128. 